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We build a "bottom-up" holographic model of charmonium by matching the essential spectral 
data. We argue that this data must include not only the masses but also the decay constants of the 
J/i/' and tp' mesons. Relative to the "soft- wall" models for light mesons, such a matching requires 
two new features in the holographic potential: an overall upward shift as well as a narrow "dip" 
near the holographic boundary. We calculate the spectral function as well as the position of the 
complex singularities (quasinormal frequencies) of the retarded correlator of the charm current at 
finite temperatures. We further extend this analysis by showing that the residues associated with 
these singularities are given by the boundary derivative of the appropriately normalized quasinormal 
mode. We find that the "melting" of the J/ip spectral peak occurs at a temperature T ~ 540 MeV, 
or 2.8 Tc, in good agreement with lattice results. 
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I. INTRODUCTION 

The fate of quarkonium states such as J/ip in the 
quark-gluon plasma has been a subject of continuous in- 
terest for both theorists and experimentalists since the 
famous proposal [l| to use the suppression of J/ ip produc- 
tion in heavy-ion collisions as a signature of the quark- 
gluon plasma formation. The screening of color forces 
should weaken or prevent cc pair binding, manifesting in 
J/ip suppression in heavy- ion collisions capable of cre- 
ating quark-gluon plasma. Making detailed predictions 
for this effect is a formidable challenge, requiring under- 
standing of such factors as cc pair production, binding, 
and survival of charmonium states in medium, as well as 
their interplay in the dynamically evolving high energy 
density medium (see, e.g., Ref. [g for review). 

A field-theoretical quantity which encodes much, but 
certainly not all, information necessary in this regard is 
the spectral function of the operator of charm current, 
jM ~ cj^c. This operator, acting on the vacuum state, 
produces charmonium 1 states such as J/ip, ip' , etc. 
By studying the response of the thermal medium to the 
action of such an operator we can learn about the fate of 
the charmonium states in this medium. The dissociation, 
or "melting" , of the charmonium states manifests itself 
in the gradual decrease with temperature of the strength 
of the response at the frequencies at which the vacuum 
state would "resonate" by creating a charmonium bound 
state. 

The current-current correlator and its spectral function 
are field-theoretically well-defined quantities, interesting 
in their own right. They can be, in principle, calculated 
using ab initio lattice Monte Carlo approach. The task of 
extracting the charmonium spectral functions from dis- 
crete Euclidean correlators measured on the lattice is dif- 
ficult, but significant progress has been made recently us- 
ing the Maximum Entropy Method (MEM) [3- 10]. The 
spectral peak corresponding to the J/ip state appears to 
be still prominent up to temperatures as high as almost 



2.5 Tc, where w 190 MeV (see, e.g, Ref. [ll| for a re- 
view) is the critical/crossover temperature at which the 
thermal QCD medium undergoes deconfining transition. 
To understand lattice results better, it is very desirable 
to have analytically controllable models which could al- 
low us to study the "melting" of the quarkonium spectral 
peak. 

The models based on a non-relativistic Schrodinger 
potential have been traditionally used for this pur- 
pose. In the pioneering work of Ref. [l3|, the effective 
temperature-dependent potential was modeled by intro- 
ducing Debye screening into the linearly confining po- 
tential. The recently studied models have used a non- 
relativistic Schrodinger potential, whose dependence on 
temperature is introduced either by input from the lat- 
tice Ts'-'^T'l , or resummation-improved perturbation the- 
ory [22.--^]. 

The non-relativistic thermal potential models typi- 
cally predict charmonium dissociation at temperatures 
approximately 1.2-1.5 Tc [lB|, which appear to be too 
low compared to the dissociation temperatures inferred 
from the lattice. Several possible resolutions of this puz- 
zle have been proposed (as reviewed, e.g., in Refs. [Tsl . 

MM)- 

Gauge-gravity holographic correspondence [27l - [29| 
provides a new field-theoretically consistent framework 
for modeling the properties of thermal strongly-coupled 
medium (see, e.g., [30,|3l| for recent reviews). Although 
ab initio calculations in QCD-proper are still not possi- 
ble within this approach, many properties of QCD, both 
in vacuum and at finite temperature can be modeled. 
"Top-down" reduction of string theory on suitably cho- 
sen backgrounds allows one to study QCD-like field the- 
ories, such as theories with flavored matter in fundamen- 
tal color representation [32l - [35| . Alternatively, reversing 
the rules of the holographic correspondence one can build 
"bottom-up" models by matching the relevant features of 
QCD, such as chiral symmetry breaking, operator prod- 
uct expansion constraints on correlators and linear con- 
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finement (sB-Sll. 

The fate of meson states in strongly coupled thermal 
medium has been a subject of many studies [H, l43l - [53| 
using a "top-down" approach to implementing massive 
flavor degrees of freedom introduced in Refs. [H, |53 |. 
Spectral functions and quasinormal modes were studied 
in Refs. [13, IH, [sO, IH, H^]- Structure of quarkonium in 
these theories was also studied in [ssj . 

It is reasonable to expect that reproducing the charmo- 
nium spectrum correctly is an important prerequisite for 
predicting its finite temperature properties. One impor- 
tant way in which the spectrum of heavy quarkonium ex- 
citations in "top-down" models differs from QCD is that 
in these models the lowest meson mass mi emerges as 
the only scale determining the masses of both the ground 
state and the excited mesons (56j : m„ ~ nrrii, n 3> 1. 
In contrast, in the heavy quark limit of QCD, the mass 
of the ground state, such as J/V'j is controlled by a pa- 
rameter (heavy quark mass) different from the parameter 
(string tension, or Aqcd) controlhng the level spacing of 
excited states: ^ ml + tiAq^j-,, Refs. [13 It 
would be interesting to find a holographic model which 
could represent these features of quarkonium spectrum 
correctly. 



II. APPROACH AND OUTLINE 

In this paper, we shall approach this problem from the 
complementary side of the "top-down" models. Taking 
the quarkonium properties at zero temperature as input, 
we shall construct the holographic dual which matches 
these properties in the spirit of the "bottom-up" ap- 
proach of Refs. [40l - [43 |. Then we shall study how such a 
model would evolve with temperature. 

The pioneering calculation of the J/ip spectral func- 
tions within such an approach has been performed in 
Ref. [H, H^l • Building on the success of the "bottom-up" 
models for light mesons [iol-li^ the authors of R ef. 
changed the scale in the "soft- wall" model of Ref. [i^l to 
match the mass of J/ijj meson, in place of the p meson. ^ 
However, such an approach suffers from the drawback al- 
ready discussed above for the "top-down" models: the 
scale that sets the level spacing of quarkonium excited 
states is the same as the scale that sets the mass of the 
ground state, in effect Aqcd ^ m-j/v the models of 
Refs. [59l - [6^ . What we would like instead is a spectrum 
which has a gap between the vacuum and the ground 
state whose scale is set by a parameter independent from 
the spacing of higher excited levels. Experimentally, the 
slope of the radial excitation trajectory for J/ijj is similar 
to that of the p mesons [13, HI] , and we want the model 



Similar "rescaling" has been used in earlier papers in order to 
calculate the charmonium dissociation temperature, Refs. [6ll . 
|62|| . but not the spectral functions. 



to match that property of QCD. 

The matching of the quarkonium mass spectrum is not 
the only new ingredient which we need to introduce. We 
point out that another important quantity to be matched 
is the decay constant of the quarkonium. Simply put, 
in order for the spectral function to be correct at high 
temperature, we should begin with a function correct at 
zero temperature, which requires matching not only the 
position, but also the strength of the J/ip resonance, i.e. 
the decay constant. 

One of the most striking qualitative consequences of 
such a more realistic model of quarkonium is a much 
more robust spectral peak of J/tp, which persists up to 
temperature 540 MeV, i.e., 2.8 Tc, in agreement with lat- 
tice studies. This is in contrast to Ref. [59] which finds 
charmonium peak "melting" at around 1.2 Tc. We take 
better agreement of our results with the lattice calcula- 
tions as evidence that the model we introduce captures 
essential features of quarkonium. 

In Section [ml to make the paper more self-contained, 
we lay out the general setup of the holographic model. 
Section |IV] discusses the features of the holographic po- 
tential we introduce to model charmonium. Section IVl 
applies the methods of holographic QCD at finite temper- 
ature to our model and, similar to Section [llll should be 
familiar to practitioners. The finite temperature spectral 
functions obtained using the proposed new holographic 
potential are presented in Section |Vll In Section I VIII 
we further our study of the thermal properties of char- 
monium by considering quasinormal modes. To facilitate 
this analysis, in Section [Villi we generalize the quantita- 
tive relation between the holographic wave-function and 
the residue at the corresponding singularity of the 2-point 
correlator (i.e., decay constant) to the case of quasinor- 
mal modes at nonzero temperature. Interestingly, this 
relationship is the same as the known one in vacuum, 
after the quasinormal mode is appropriately normalized. 
We derive the expression for such a norm in Section [Villi 
In Section IIX[ we use the quasinormal modes and the 
residues to analyze the charmonium spectral function 
more quantitatively. We conclude in Section |X] with a 
summary and a discussion of the results. 



III. SETUP 

We shall focus on the two-point correlation function of 
the heavy quark (charm) current, = cj^c, and define 
the retarded Green's function: 

GR{q) = -i Jd''xe''^-e{xo) {[J^{x),J^iO)]), (1) 

where is a component of the current orthogonal to the 
4- vector q^: q^Jf^ — and (. . .) denotes thermal average. 
For simplicity, we shall consider only the case q = 0, i.e., 
q = (a;,0). Thus will denote any spatial component 
of J^, e.g., J^. We shall also define Gr{uj) = Gr{uj,0). 
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The spectral function is defined, for real values of uj, as 



p{lo) = p{uj, 0) = - Im Gr{uj). 



(2) 



Among standard properties of jo(a;), which can be shown 
using spectral representation, is that /o(a;, q) equals the 
Fourier transform of {[J^{x), J-'"(0)]), and that p{uj)/u! > 
for all u!. 

In the spirit of the holographic approach, we shall as- 
sume that the generating functional of the heavy quark 
vector current can be represented by the effective ac- 
tion obtained by integrating over a bulk 5D gauge field 
Vm (dual to the current) at given fixed boundary values 
(equal to the source of the current). The action for the 
5D gauge field is given by 



S 



(3) 



where is the 5D gauge coupling and Vmn — QmVn — 
dNVM- The 2-point current correlator is given by the 
linear response of the field Vm to an infinitesimal pertur- 
bation of its boundary condition. 

We choose the conformally flat representation for the 
5D background metric guN with 4D Lorentz-isometry: 



gMN dx^^dx^ = e2^(^) [rf^^dx^'dx" - dz^] , (4) 



where rj^i, = diag(l, — 1, — 1, — 1) is the Minkowski met- 
ric tensor. The effect of confinement is represented by 
the non-trivial background profile of the scalar field $ 
in Eq. ([3]) in the same way as it is done in the soft- 



wall model with dilaton background in Ref. [42|. We 
shall assume that the physics associated with the mass 
of the quark and chiral condensate operator cc can be also 
represented as a contribution to the background field 
Such a contribution is essential if the model is to be ex- 
tended to describe different flavors, which must have dif- 
ferent holographic backgrounds, depending on their re- 
spective masses. In a putative top-down approach this 
contribution could arise through the interaction of the 
gauge held with the scalar ("tachyon") held dual to the 
quark mass operator cc, e.g., by a mechanism explored 
in Refs. 

Following the rules of the holographic correspondence, 
we shall calculate the generating functional for correla- 
tion functions of the heavy-quark current by evaluating 
the action at its extremum for given boundary conditions. 
The extremum is given by the solution of the equations 
of motion, which in V5 = gauge read: 



F 



(5) 



where V is any of the three components V±^ of V^(g, z) 
transverse to 4-vector g'' [q^ — 0) and 



(6) 



Discrete values of 



\, for which the Eq. ([5]) pos- 



sesses a normalizable solution V = Vn{z) satisfying the 



boundary condition V|z=o = 0, correspond to the masses 
nin of the charmonium states, n = 1, 2, . . . = J/ip, ip' , .... 
We normalize such solutions as 



dze^'^'^Vniz)^ = 1 



(7) 



Back on the held theory side, the decay constants /„ 
are defined via matrix elements of the charmonium cur- 
rent between vacuum and a given vector charmonium 
state n with polarization as 



\J^{0)\n) = /„m„e^ 



(8) 



Holographic correspondence relates these constants to 
the (second) derivative of the normal mode Vn (see 
Refs. 0113): 



fn 



v'„{z)e 



B{z) 



(9) 



As one can see from Eqs. ©-(H]), the masses and 
decay constants of the charmonium states are deter- 
mined by the combination B of the gravity {A) and dila- 
ton/tachyon ($) backgrounds ([6]). In the spirit of the 
bottom-up approach, we shall choose the function B so 
as to satisfy a number of phenomenological and field- 
theoretical QCD constraints. We assume that such a 
background arises dynamically, but do not attempt to 
model the corresponding dynamics. 

We require that the ultraviolet behavior of the 
current-current correlator is conformal, viz. Gji{uj) 

log(— w^) as — )• —00, which translates into: 



^0 _i 
— > z 



and matches that of QCD, which fixes [371 l4i 
gl = UiryN, . 
By performing a Liouville transformation. 



(10) 



(11) 



(12) 



we can bring Eq. ([5]) to the canonical Schrodinger-like 
form 



- d'^'i'/dz^ + U{z)^ = qH 
with the holographic potential given by 



(13) 



(14) 



We shall now choose B or, equivalently U, to match 
not only the mass of J/ip, but also the mass of tfi' as well 
as (and more importantly) the decay constants of J/^p 
and ip' . 
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IV. HOLOGRAPHIC POTENTIAL 

In order to explain the role and motivation for the 
features of the holographic Schrodinger potential which 
we are led to introduce in this paper, we shall add those 
features one at a time. 

To begin with, it is easy to implement the two-scale ex- 
citation spectrum, correcting the drawback of Ref. [s^, 
by adding a constant, c^, to the potential. Taking the 
potential from the soft- wall model, C/(a): which repro- 
duces the equidistant mass-squared spectrum controlled 
by parameter a, we thus consider 



f/(a,c) - U(a) + - 3/(4z2) + [a" zf + 



(15) 



For this potential, the mass of the ground state is con- 
trolled by c, independent from a? 

The shifted potential in Eq. (fTS]). however, has an im- 
portant drawback, which is significant as far as the finite 
temperature properties of the quarkonium are concerned. 
Since those properties are encoded in the spectral func- 
tion of quarkonium at finite temperature, it is clear that 
to get those properties right on a quantitative level, one 
must begin with correct spectral function at zero temper- 
ature. The spectral function at zero temperature consists 
of a series of peaks, representing the quarkonium states. 
These peaks are characterized not only by their position 
(masses) but also by their strengths, i.e., the decay con- 
stants of the quarkonium states. If we take the model of 
Ref. [13, we find that the decay constant of J/ip is un- 
derpredicted by 20% and that of tp' is overpredicted by 
15% -- see Table m For the shifted potential in Eq. ([T5l). 
the situation is much worse: the decay constants are un- 
derpredicted by a factor of order 2 — 3. 





Experiment 




Observable 


(MeV) 


(MeV) (MeV) 




3096 


3096* 3096* 




3685 


4378 3685* 


fj/ip 


416 


348 145 


ftp' 


296 


348 173 



TABLE I: Comparison of the masses and decay constants 
obtained using the scaled soft-wall (as in Ref. '59']) and 
shifted soft-wall potentials in Eq. (|15p with experimental val- 
ues, Ref. [53. Analytic formulas from Ref. [12] are used: 



m„ = 4a n -I- c and /„ — V8na /{gsnrin), n = 1,2. The 
observables which are fitted to determine parameter a in U(^a.) 
and a, c in U(a.c) are marked by asterisks. Both models signif- 
icantly underpredict the decay constant of J/t/'- The shifted 
potential f/(a,c) is worse because the parameter a needed for 
the fit is smaller. 



What should be done to the potential ([T5t to increase 
the value of the J /if; decay constant? To understand the 
solution to this problem, one should recall the relation- 
ship between the decay constant and the rate of change 
of the holographic wave- function at the boundary z — >■ 0, 
Eq. da]): 



fn 



■<(o) 



(V^^n)' 



(16) 



where we used Eq. (1101) and the normalizable solution ■(/;„ 
to Eq. ((T3)) . = e^^^Vn- That means the steeper the 
wave- function tpn is near the boundary z = 0, the larger 
is the decay constant. For example, increasing parame- 
ter a would make the function steeper, by squeezing its 
support. But that will unacceptably increase the radial 
trajectory slope, drn^/dn, which we want to preserve. A 
simple and natural solution is to provide an "incentive" 
for the wave function to be larger only in a narrow re- 
gion close to the boundary, thus forcing it to approach 
the boundary value V'n(O) = steeper. This can be eas- 
ily achieved by creating a "dip" in the potential at small 

Z ~ Zd- 

What should be the shape of the dip at zjl In this 
paper we shall consider a minimalistic example of the po- 
tential necessary to match the quarkonium spectral data. 
We shall use the simplest approximation for a narrow dip: 
a negative Dirac delta function.'^ 

Rather than superimposing the delta function onto the 
potential C/(a.c) in Eq. p3)l . we first observe that, since 
the potential undergoes a qualitative change in behavior 
at z = Zrf, manifested in the delta function, there is no 
reason to expect that it is described by a function contin- 
uous across Zd- A more natural, and convenient, choice 
is to consider the following piecewise analytic function: 



U{z) = — 6{zd -z) + {{a'zf + c^) e{z - Zd) 

-a5{z-Zd), (17) 

which represents the necessary large and small z behav- 
ior. There are 4 parameters in this potential, and we use 
them to fit 4 experimental data points: the masses and 
the decay constants of J/ij} and ij}' in Table ID We find 
(see also Fig. [T|) 



a = 0.970 GeV, c = 2.781 GeV, 
a = 1.876 GeV, Zd = 2.211 GeV. 



(18) 



One can see the effect of the delta function on the holo- 
graphic wave functions explicitly as shown in Fig. [5] As 
expected, the wave functions, varying relatively smoothly 
for z > Zd, "dive" steeply towards their boundary value 
i'niO) — once z < Zd- 



^ The first term is the potential for a conformally invariant theory. 
It can be, curiously, viewed as a centrifugal term for a radial 
Schrodinger equation in 2d with the radial quantum number m = 
1 (coinciding with the spin of the mesons we study [i^l ). 



Replacing the delta-function with a narrow square well produces 
similar results. 
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U(z) (GeV^) 




FIG. 1: Holographic Schrodinger potential at T = given 
by Eqs. (|17p and (|18[l . The negative delta-function a,t z — 
Zd is indicated by an arrow. Dashed horizontal lines show 
the two lowest eigenvalues, m\ and ml {J/i^ and ip'). The 
discontinuity &t z = Zd plays less significant role compared to 
the delta-function itself. 



iA(z) 




-0.5 - s 
-1.0 - 

FIG. 2; The first two holographic wave functions il'i and i/)2 
(corresponding to J/?/; and ijj' mesons) for the potential in 
Fig. [T] The kink &t z = Zd is due to the dip (delta function) 
in the potential. 



The potential in Eq. (|17l) can be, certainly, improved 
by applying further constraints, e.g., fitting the masses 
and decay constants of the higher excited states, and also 
by applying constraints which follow from the operator 
product expansion of the heavy quark current-current 
correlator [G^ beyond the leading order."* We shall de- 
fer such improvements to future work, and concentrate 
here on the consequences of the two major properties of 
the potential we introduced: the shift, providing the two- 
scale level pattern, and the dip at small z, supporting the 
necessary size of the decay constants. 

From the point of view of the nonrelativistic model of 
quarkonium, the effect of increasing the decay constant 
can be understood as the effect of decreasing the size 
of the bound state: the decay constant is proportional 



* In fact, our removal of the term from the region of small z, 
achieved by using piecewise potential l|17| l. is motivated, in part, 
by such considerations. 



to the probability of finding quark and antiquark at the 
same point, which is larger for a more compact state. It 
is reasonable to expect that a more compact bound state 
survives up to higher temperatures. This is indeed the 
result we find (see Section rvT| . 



V. FINITE TEMPERATURE 

Once the zero-temperature parameters of the model 
are fixed, we wish to determine the spectral function of 
the charm current at finite temperature, as defined on 
the field theory side by Eqs. ([J) and 

On the holographic side, finite temperature correla- 
tion functions of heavy quark current are determined by 
a similar action to Eq. but on a different, temper- 
ature dependent, background $ and qmn, whose main 
feature is the presence of a black-hole horizon at some 
value oi z = Zh- In a "top-down" approach that back- 
ground would itself be an extremum of the same grav- 
ity /dilaton/tachyon action. We shall not attempt to 
model the dynamics of these background fields here, leav- 
ing it to further work. Lacking the equations of motion, 
which would determine the background for each temper- 
ature, we shall instead make a minimalistic assumption 
about the behavior of this background. 

We choose the following representation for the metric 
consistent with the 3d Euclidean spatial isometry, corre- 
sponding to finite temperature: 

ds" ^e^^^'^[hdt' -dx" -h-^dz''] . (19) 

If the function h{z) has a simple zero, viz. h{zh) — 0, 
the space described by possesses an event horizon 
at z = z/j. The temperature T corresponding to this 
background is related to z^ as 

T=^\h'{zn)\ . (20) 

We assume the simplest ansatz for h, interpolating be- 
tween h{Q) — 1 and h{zh) — with the power z^ dic- 
tated by the dimension of the operator of the energy- 
momentum: 

h=l- {z/zhf , (21) 

which turns out to be the same as in the familiar AdS^ 
black hole solution to Einstein equations with negative 
cosmological constant. The temperature T correspond- 
ing to this background is related to z/j as 

Zh = (^T)-i . (22) 

As already noted, in principle, we should expect the func- 
tion A{z) as well as the background <I>(z) to depend on 
the temperature. Since we do not attempt to model the 
background dynamically, we have little choice but to ne- 
glect the dependence of the background B in Eq. ([5]) on 
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T. Our purpose here is to describe the effect of the tem- 
perature semiquantitatively. The presence of the black- 
hole horizon in the metric (|19p is the main property that 
our ansatz is aimed to capture. 

Our approach to introducing temperature dependence 
here is similar to that of Ref . [59| . ^ In addition to making 
our model simpler, this will allow us to see the effect of 
the features of the potential we introduce (specifically, 
the "dip" ) , in comparison with Ref. [S^ , more clearly. 

Therefore, the equation we solve to determine the cor- 
relator of the heavy quark current at finite temperature 
is given by 



d,{he''d,V)+u^h-'e''V = 



(23) 



where function B{z) is determined by solving Eq. for 
given U{z). The retarded correlation function is given, 
according to the well-known prescription of Ref. [68., .69,] , 
by 



GR{oj)^~\he''V\z,uj) 

95 



1 V'{e,u;) 



9l 



(24) 



where e — is an ultraviolet regulator and T^(z, w) is the 
solution of Eq. ([23| with boundary conditions: 



(25) 



y(w,e) = l ; 

V{lo,z) C{u:){l - z/z;,)-^"/(4-^) , 



where, for each w, C{uj) is a free constant determined, 
if needed, by solving equation (|23l) . The spectral func- 
tion can be calculated as the imaginary part of Gr{uj) 
given by Eq. (p4)) . Alternatively, one can observe that 
lm{he^ V* dzV) is z-independent and evaluate it at z = 
Zfi, instead of z = e: 



(26) 



which demonstrates explicitly that p{uj)/uj > condition 
is a built-in property of the model. 

The Green's function obtained using Eq. is loga- 
rithmically divergent as e — >■ 0. This logarithmic depen- 
dence on the UV regulator e matches the expectation on 
the field theory side. It does not concern us here because 
the imaginary part of Eq. (p4)) . i.e., the spectral func- 
tion given by Eq. (P^ . is finite, as it should be in the 
corresponding field theory. 

The solutions of Eq. (^5]) are easier to study by apply- 
ing the Liouville transformation 

*(C(z)) - e^(^)/2 V{z); C{z) - dz'/h{z') , (27) 



which yields the equivalent Schrodinger equation for the 
function 4'(C): 



(28) 



^ Similar approach was used earlier in Refs. [69 . \6l\ to study light 
flavor mesons at finite temperature. 



where the potential is given by 



B"{z) ^ fB'{z)Y B'{z)h'{z) 
^ 2h{z) 



h{zf . (29) 



At T = this equation becomes Eq. plj) with C, — z. 
We use Eq. ([H]) with U{z) given by Eq. ^} to find 
B{z), which then, upon substitution into Eq. p9t gives 
us C/t(C(^))- 



VI. RESULTS 

For the purpose of presenting the results, it is useful 
to observe that in a theory with conformal behavior at 
short distances, one can expect that 



p(w)/a;2 



dimensionless const. 



(30) 



This large uj behavior is dominated by the contribution 
of the unit operator in the OPE of the product of the 
currents. In QCD, the dimensionless constant can be cal- 
culated, due to the asymptotic freedom, p(cLi)/a;^ ~ — 



Nc/{24:Tt), and it is matched in the holographic model by 
Eqs. (Uni), (HH): 



p{uj) 



25i 



(31) 



Equation motivates the use of the rescaled spectral 
function: 



P/^^ 2g| p 



(32) 



such that p(oo) = 1. The spectral functions of the charm 
quark current obtained in our model at three representa- 
tive temperatures are shown in Fig. |3l 

We observe that a pronounced J/ V' spectral peak sur- 
vives at least beyond T = 400 MeV « 2Tc. This is in 
agreement with lattice results discussed in the introduc- 
tion. In comparison, the model in Ref. (59| , predicts the 
disappearance of the peak already above T ~ 1.2 Tc. We 
can also see that the J/ip peak has mostly disappeared 
at T = 600 MeV. 

Another interesting observation is that the peak at zero 
UJ grows with temperatures. This peak is related to DC 
conductivity a by Kubo formula 



a — lim p{lj)/ijJ 



(33) 



It is easy to show, using Eq. (pS)) and the fact that V =1 
is the solution to Eq. ([231) at w = 0, that (c.f. Ref. [70]) 



(34) 
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FIG. 3: Rescaled spectral function, Eq. at T = 200 (solid 
curve), 400 (dashed curve), and 600 MeV (dotted curve) in 
the holographic model of this paper. 
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FIG. 4: Holographic potential in Eq. ^ at T = 200 and 400 
MeV. See also Fig. [l] 



On the field theory side, the growth of the charm con- 
ductivity is a natural consequence of the dissociation of 
quarkoniuni bound states. 

It is also instructive to look at the temperature evolu- 
tion of the holographic potential. Comparing the poten- 
tials at r = 200 and 400 MeV in Fig. HI one can see that 
the barrier protecting the (quasi-)bound state from de- 
cay becomes thinner, and the depth of the delta function 
dip becomes smaller with increasing temperature. 

In order to obtain a more objective and quantitative 
description of the evolution of the size and shape of the 
spectral peak, we shall now take a look at the quasinor- 
mal modes defined by equation (1231) with outgoing wave 
boundary condition from Eq. (|25p . 



VII. QUASINORMAL MODES 

In order to study the properties of the peak in the 
spectral function, corresponding to J/ip, we shall use its 
relation to one of the poles of the function Gii{u!) in the 
complex plane. Such poles correspond to solutions V = 
Vn of the equation (|23p satisfying boundary conditions 



(35) 



similar to the b.c. in Eq. (1251) at the horizon z — z^, but 
not at the boundary z = e. Unlike the b.c. in Eq. (|25l) . 
the b.c. in Eq. (I55|) are homogeneous and specify the 
solution only up to an overall normalization constant. 
We describe the most suitable normalization condition 
for Vn in Section IVIIII 

Such solutions u„ are called quasinormal modes due to 
their similarity to the normal modes of the T = equa- 
tions. Similarly to normal modes, the quasinormal modes 
exist only for a discrete set of frequencies a;„. However, 
due to the nature of the b.c. in Eq. ([55]) . these frequencies 
are complex. 

Quasinormal modes in the context of gauge/gra vity 
duality have been studied extensively, see, e.g., Refs. O- 
pTSj . and in particular, in application to the fate of mesons 
at finite temperature in Refs. [13, [H, [13, [H, [13]. We 
follow a standard numerical method to obtain the val- 
ues of the quasinormal frequencies as well as the corre- 
sponding eigenfunctions of the mode equation. In sum- 
mary, the main challenge arises due to the fact that, 
since Ima;„ < 0, the boundary condition selects the so- 
lution dominant as z ^ Zh- Thus, numerical integra- 
tion of the equation must be carried out with precision 
much better than (1 — z/z;i)~^™""/(^'^-^\ to make sure 
the "wrong" solution is not mixed in. This effectively 
prevents the numerical integration from reaching z — Zh- 
The standard solution to this problem [tI], and the one 
we adopt, is to build the solution in a small but finite 
interval [z/i(l — S),Zh] using a truncated Taylor (Frobe- 
nius) expansion. The numerical integration can then be 
carried out starting from the point Zh(l — 5) away from 

Zh- 

The trajectories of the first two quasinormal frequen- 
cies on the complex plane are shown in Fig. [51 One can 
see that the imaginary part of the second mode W2 in- 
creases much faster with temperature than that of wi. 
This corresponds to the observation that the ip' peak 
"melts" very early, in accordance with the sequential dis- 
sociation scenario [77| . 

We can demonstrate the relationship between the 
quasinormal modes and the peak more clearly by cal- 
culating the residue r„ of Gr in the pole corresponding 
to the mode n 

r„ = liia {uj - ujn)Gii{uj) . (36) 

a;— >-cJti 

It is easy to see that at T = 0, the residue is related to 
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FIG. 5: Quasinormal modes uji (J/tp) and aj2 {tp') for a se- 
quence of temperatures: T = 100, 200, 300, 400, 500, 600 MeV 
for cji and T = 100,200,300 MeV for L02. Colors (online) 
correspond to temperatures. 



the r = decay constant in Eq. © as 

Jn"''n 



x2 

2w„ 



(r = o) 



(37) 



We shall define the contribution of the pole to the 
rescaled spectral function (15^ as 



Pn{iL>) = Im 



where 



TT LUf, 



(38) 



Due to the well-known property of the Green's func- 
tion, Gr{lu)* = Gr^—uj*), the poles must come in pairs, 
w„ and uj-n = ^'^nj with residues r„ and r_„ = — r* . 

In Section IVIIIl we show how the quasinormal mode 
could in fact be normalized and derive a useful expression 
for the residue r„ in terms of the corresponding quasi- 
normal mode. The result is an extension of the T = 
equation 

To demonstrate the usefulness of the expression for r„, 
given by Eqs. ([5^ . (IMl) we derive below, we show in Fig.[S] 
the spectral function p, together the contribution of the 
pole and its "mirror," pi + p^i, as well as the result of 
the subtraction p — pi — p_i which, as expected, leaves 
only monotonic background around the location of the 
peak, uj ~ Reoji. ^ 



VIII. NORMALIZATION OF THE MODES AND 
THE RESIDUES 

In this Section, we describe how the residue r„ for 
a pole at quasinormal frequency aj„ can be calculated 



Although, in the case we consider, the contribution p_i of the 
"mirror" pole at — is almost negligible in the vicinity of the 
peak at til Ri RetJn, we keep this contribution just to preserve 
the symmetry of the spectral function 




(J (GeV) 



FIG. 6: The scaled spectral function at T = 200 MeV. The 
dashed line is the contribution of the first quasinormal mode, 
pi + The dotted line is the result of subtraction of that 
contribution from the spectral function. 



given only the corresponding quasinormal mode Vn{z). 
We shall show that the expression is similar to the zero 
temperature result, Eq. for the decay constants /„. 
In order to do that, we would need to be able to gener- 
alize the normalization condition Eq. ([7]) at T = to the 
case of quasinormal modes at finite T. We shall indeed 
obtain this generalization as a byproduct of our deriva- 
tion of r„. 

Similar to the derivation of the T = result 
in Ref. [i^, we shall introduce the Green's function 
G{ijj; z, z') for the mode equation 



(he^d.G) + uj^h-^e^G = S{z - z') 



(39) 



with (homogeneous) boundary conditions similar to 

Eq. (I2SD 



(40) 



G(a;;e,z') = ; 

G(c.; z, z') C(c., z'){\ - z/z,)— 



where, C(a;, z') is an arbitrary function of uj and z'. It 
can be easily checked that 



y(c^,z') =-z"i9,G(t^;z,z')U=. 

satisfies defining equations and (I^SI) . 

The Green's function G can be constructed as 



G(aj; z, z') 



u(w, z<) i;(a;, z>) 



(41) 



(42) 



where Z{< >} = {min, max}(z, z'); u and v are two solu- 
tions of equation ((23|) each obeying separately one of the 
two boundary conditions in Eq. (|40|) : 



u(uj, e) = ; (43a) 
t;(w, z) - (1 - z/z/O"'"/^^''^^ for z^Zh . (43b) 



Both boundary conditions (|43p are homogeneous and 
thus define u and v only up to overall normalization con- 
stants. However, the normalization of either u or w is 
canceled out from Eq. due to the Wronskian 



[u, v\=u dzV — V dzU , 



(44) 
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in the denominator. Note also that the denominator in 
Eq. (021) is independent of z: dz(he^\u,v\) = 0. 

The zeros of the Wronskian [m, v\ as a function of uj are 
the poles of the Green's function (P^ . The corresponding 
complex value of w is the quasinormal frequency, w„ . The 
residue associated with such a pole is given by 



lim (w — a;„)G(a;; z, z') 



z<) ti(w, z>) 



(45) 

Since he^di^[u, v] is independent of z, we can calculate 
it by setting z to a convenient value, which happens to 
be z = e. At that point e) = and di^u{uj, e) — 0. 
Furthermore, since at w = aj„ the function w is a constant 
multiple of u, also v{ujm e) = 0. This means. 



(at uj = ujn). (46) 



Differentiating Eq. (f23|) obeyed by v with respect to ui, 
we obtain the equation for d^iV 

dz {he^dz{duv))+Lu^h-^e^{du:v)+2ujh-^e^v = 0. (47) 

We multiply Eq. (gT]) by m, then take Eq. (^31) obeyed 
by u and multiply it by d^v, and then subtract the two 
resulting equations to find 



-2u}h ^e^uv 



(48) 



Integrating Eq. P5|) over z from z = e to z = z;j(l — (5) 
and using Eq. we find (at u = w„) 



he^d.^Va, v] = 2uj 



■ e uv 



(49) 



z=Zh(l-(5) 



As ^ — >■ 0, we can (provided — Imaj„ < 2tiT, see below) 
replace v in the last term by its asymptotics (j43bp . Thus 



OuV > 



4ttT 



log(l - z/zh)v, 



(50) 



and therefore 



z={l-5)zh 



^ — [u,v\ + —uv . (51) 



At a; = w„, u ^ V and [m, u] = 0. Thus we can write, 
for 6^0, Eq. (gH) as 

ae 0^^1*1^ =261; / — e uv + ie uv\ . 

I- ' J J ii i2=(i-i5)2h 

(52) 

If we choose the normalization of u and v aX u — uin in 
such a way that u ^ v ~ Vn and the r.h.s. of Eq. (j52p 
equals 2a;„, we find from Eqs. (gS]), dH]), ^ and ((M)) : 



1 «(^)/e)^ 
gl 2ujn 



(53) 



where a;„ is the quasinormal frequency and w„ is the cor- 
responding quasinormal mode normalized as 



lim 

5^0 



I 

2u)r. 



z = {l-S)zh 



= 1 . 



Note that, at T = 0, the normalization condition in 
Eq. (|54l) reduces to the normalization in Eq. ([7|) for the 
normal mode. The expression for the residue in Eq. (j53p 
also reduces to dH) via Eq. ([37)) . 

It is important to realize, however, that is involved 
in Eq. ([54| . not \v'^\, and that both terms in Eq. ([54)) 
are complex. In particular, that means that complex 
condition in Eq. (IMl) fixes both the magnitude and the 
phase of u„. 

One could also note that J{dz/h)e^v^ = f d^-ip^, 
where C and -0 are the Schrodinger variables defined in 
Eq. ([?7)) . Written in terms of C, and ■0, the normalization 
condition in Eq. ([7]) is equivalent to the one introduced 
by Zeldovich in the study of the quasi-discrete levels in 
quantum-mechanical decays [tI], and used in the theory 
of quasinormal modes of cavities [79| . 

Equation ()54p has a drawback, which becomes im- 
portant when — Ima;„/r is sufficiently large. Since 
Imw„ <0,Vn-^{l- z/zft)-'"/(4-^) — ^ 00 as z — z/j. 
The last term in Eq. (|54p serves to cancel the diver- 
gence in the integral. However, the remainder is fi- 
nite only if — Inicjn < 2ttT. This follows from the 
Frobenius expansion of the solution v near z — z^' 
viuj, (1 - 5)zh) ^ ,S-^'^/(4'^^)(l + 0{S)). Substituting into 
Eq. ([54)) . we find that after cancellation of the diver- 
gent terms J^^^/C^tt^)^ remaining terms are of order 
^-zu/i27TT)+i g^^^ vanish as S for -Imw„ < 2TrT. 
However, for — Imw„ > 2ttT, the remaining terms are 
still divergent as (5 — > 0. This drawback can be cured 
by returning to Eq. (|49p and keeping more terms in the 
asymptotics of dc^v in Eq. (150)) . In practice, we do use 
Eq. ()49p . and calculate the last term on the r.h.s. us- 
ing the truncated Frobenius expansion (which we have 
to generate anyway in order to obtain the solution of the 
mode equation near z/i, as explained in Section [VIip . 



IX. RESULTS OF QUASINORMAL MODE 
ANALYSIS 



The analysis of quasinormal modes allows us to study 
the fate of the J/ip more quantitatively. Rather than 
fitting the spectral function p, as done in Ref. [59j , we 
can read off the peak's height and width directly from 
the values of the residue and the imaginary part of the 
quasinormal frequency obtained by solving Eq. (I23p . with 
b.c. ([55)) and using Eq. (155)) and ([5^ . We define the peak 
height Hn as 



CJ — Rc UJn 



Ima;„ 



(55) 
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FIG. 7: Temperature dependence of the normalized height 
of the J /tp peak Hi defined in Eq. (|55|) in the holographic 
model. For comparison, the same result using the potential 
of Ref. [59l |. i.e., U(a) of Eq. (|15|l . is shown as a dashed line. 
The dotted horizontal line shows asymptotic value of p„. 

where p„(w) and f„ are given by Eq. ([38| . The width of 
the peak can be defined, as usual, as 

r„ = - Imw„ . (56) 

It is instructive to compare our results with those ob- 
tained by applying the same quasinormal mode analysis 
to the scaled soft-wall model [595. 

The dependence of the J /i) peak height Hi on T is 
shown in Fig. [T] Using Hi, one can define a useful cri- 
terion for identifying the "melting" temperature, which 
quantifies reasonably well the qualitative perception. We 
can define the disappearance of the peak as a point at 
which the height of the peak as defined by Eq. ([55)1 be- 
comes smaller than the asymptotic value of the back- 
ground, i.e., p(oo) = 1. 

According to this criterion, the peak survives until 
temperatures about 540 MeV, i.e., about 2.8 Tc. In con- 
trast, in the rescaled soft- wall model of Ref. [59| the peak 
height drops quickly and disappears at about 230 MeV, 
or 1.2 Tc- The peak is more robust in our model because 
of the delta function in the holographic potential U in 
Eq. (|17l) . which can additionally be checked by reducing 
a to 0. 

The dependence of the J /^jj peak width Ti (in units of 
2tiT) is shown in Fig. [5] Interestingly, the peak begins 
to broaden somewhat earlier in our model (T ^100 MeV 
vs T ~150 MeV). This can be understood as the conse- 
quence of the fact that our potential has a "softer wall" 
at large z, and thus is sensitive to lower temperatures. 
However, even after the soft wall is "melted" at about 
200 MeV, the dip at small Zd continues to support the 
J/ij} peak. On the contrary, in the rescaled soft model 
Ref. [5!| , once the wall has "melted" , no feature remains 
in the potential to support the narrow quasi-discrete J / ip 
state. 

It is also interesting to note that for T — ^ oo the width 
ri/(27rT) 1, which corresponds to the value of the 
n = 1 quasinormal frequency, uj„ = 27rTn(±l — i) (see. 



0.6 - ^ - ' / 

0.4 - ^ ' ' ^^■'^^^ 

/ ^ — 

0.2 ^-y 
/ 

V 

0.0 I T(GeV) 

0.2 0.3 0.4 0.5 0.6 0.7 

FIG. 8: Solid line: the width of the J/i/' peak in units of 
27rT as defined by Eq. (|56(l . For comparison, the dashed line 
shows the result using the potential from Ref. [s^, i.e., V^a) 
in Eq. (fTS]) . 

e.g., Ref. [tI]) in the AdS black- hole background. This 
agrees with the expectation that at asymptotically large 
temperatures the theory should appear almost conformal. 



X. SUMMARY, DISCUSSION AND OUTLOOK 

In this paper, we attempted to build a holographic 
model of charmonium using the "bottom-up" approach. 
As was done in Ref. , we required the matching of the 
UV behavior of the charm current correlator to that of 
QCD, which is known due to asymptotic freedom. This 
fixed the small z behavior of the holographic potential. 
Similarly, the large z behavior was fixed using the condi- 
tion that the meson spectrum is asymptotically equidis- 
tant in mass squared, as in the soft wall model of Ref. (i^. 

We then observed that, in QCD, the asymptotic spac- 
ing of levels is controlled by a parameter, Aqcd, theo- 
retically independent from the mass of the ground state, 
fn.j/tl)- We incorporated this feature by a shift in the 
soft wall model holographic potential. At this point our 
model became different from that of Ref. [s^ , where both 
scales were controlled by one parameter. 

We then required that not only the spectrum of the 
meson masses is matched to QCD phenomenology, but 
also the values of the decay constants. This requirement 
led us to introduce a "dip" in the holographic potential. 
We modeled this feature minimalistically, by a delta func- 
tion. 

A possible way one could anticipate such a feature from 
a putative "top-down" construction is by observing that 
the existing models describing heavy mesons do indeed 
have potentials with the characteristic scale of the width 
of the well (flavor brane extension in the z coordinate) 
of order of the inverse meson mass Zd ^ m^^ . From 
that point of view, the potential we need differs from 
these models by a slower rising "soft wall" at large z, per- 
mitting the wave-function to "leak" significantly beyond 
scales of order inverse meson mass m^^, and reaching 
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up to the confinement scale, such as confinement radius, 
-Rconf ~ ^QCD ' illustrated in Fig. [HI 




FIG. 9: A sketch emphasizing the features of the holographic 
Schrodinger potentials. Left - the "top-down" model poten- 
tial (see, e.g., Eq. (3.2) in Ref. [5l[) has one scale, m, char- 
acterizing both the ground state mass and the level spacing. 
Right - a potential with separate scales for the mass gap, m, 
and for the level spacing, Aqcd- Due to slower rising wall 
of the potential at large z, the holographic wave function is 
allowed to "leak" beyond the scale of . In this paper the 
"dip" at small Zd ~ is approximated by a negative delta 
function. 

It is easy to see that such a modification of the po- 
tential is what accounts for the broadening of the J/ip 
peak below the dissociation temperature in our model. 
In contrast, the "top-down" models would predict no 
broadening below dissociation temperature, though this 
can change if instanton-like tunneling effects are consid- 
ered fil. 

Having thus constructed a potential with 4 free pa- 
rameters, we fit those parameters to the 2 masses and 2 
decay constants of J/^p and ip'. 

Obviously, the "ascetic" potential which we introduced 
can be further improved by matching more than just the 
2 states in the charmonium spectrum. The problem of 
reconstructing the potential from the spectral data is a 
standard inverse problem. The fact that the spectrum 
together with the derivatives of the eigenfunctions (ini- 
tial velocities, or normalizing constants) uniquely spec- 
ifies the potential for a Dirichlet problem [81|, suggests 
that matching decay constants is a natural part of such 
"bottom-up" approach. 

An interesting related observation is that the same 
spectral data can be obtained, in principle, from a non- 
relativistic model of quarkonium. It would be interest- 
ing to know if a more direct relationship exists between 
the holographic potential and the potential in the non- 
relativistic Schrodinger description of quarkonia in the 
heavy quark limit. 

It might be also interesting to consider the constraints 
on the potential arising from the operator product ex- 
pansion of the current correlator [65|- By quark-hadron 
duality these constraints can substitute the spectral data. 



It would be interesting to explore the finite tempera- 
ture consequences of these constraints and compare to 
the studies of Ref. 

Having specified the model at zero temperature, we 
consider the evolution of the spectral function of the 
charm current with temperature. Here too, we make a 
minimalistic assumption about the temperature depen- 
dence of the background, similar to the previous study of 
Ref. (59| . The necessity for the assumption comes from 
the fact that we do not model the dynamics of the back- 
ground itself. This can be remedied in a more compre- 
hensive model which we leave to further work. We believe 
that the minimalistic assumption about z-dependence of 
the temperature factor h and the scale factor is suffi- 
cient to make semiquantitative predictions at the level of 
precision generally expected from such a model approach. 

One of the new tools we introduce to study properties 
of the finite temperature Green's function is based on the 
relationship between the residue in the pole and the holo- 
graphic wave function, which we derive, Eqs. (j53p . (j54p . 
This relationship is a nontrivial generalization of the zero 
temperature expression for the decay constants ^40j . This 
allows us to study not only the position of the pole in the 
complex plane, corresponding to the quasinormal mode, 
but also the strength of the contribution of this pole to 
the spectral function. In particular, this allows us to de- 
vise a simple quantitative criterion for the disappearance 
of the peak based on its height relative to the asymp- 
totic level of the background. Since the dissociation of 
quarkonium and disappearance of the peak is not a sharp 
transition, a criterion of this type might find its use in 
comparing results of different approaches in a more uni- 
form and objective way. 

Finally, we find that the new features of the holo- 
graphic potential, such as the "dip" , which we introduce 
to model quarkonium more realistically, strengthen the 
robustness of the J/ V' peak, allowing it to persist out to 
temperatures of about 540 MeV, i.e., 2.8 Tc, according 
to our criterion. This is in good agreement with lattice 
studies, but not with most non-relativistic Schrodinger 
models of the quarkonium or the holographic model of 
Ref. [1^, which typically predict dissociation at about 
1.2 Tc- The comparison with the lattice results can be 
made even more direct by calculating the Euclidean cor- 
relators of the current in the holographic model. We defer 
this to future work. 
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